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METHOD AND APPARATUS FOR POWERED 
DESCENT GUIDANCE 

CROSS-REFERENCE TO RELATED 

APPLICATIONS 5 

This application claims benefit of the following U.S. Pro- 
visional Patent Applications, each of which is incorporated by 
reference herein: 

U.S. Provisional Patent Application No. 61/201,866, 10 
entitled “MINIMUM LANDING ERROR POWERED 
DESCENT GUIDANCE FOR PLANETARY MISSIONS,” 
by James C. Blackmore and Behcet Acikmese, filed Dec. 16, 
2008; 15 

U.S. Provisional Patent Application No. 61/164,338, 
entitled “MINIMUM LANDING ERROR POWERED 
DESCENT GUIDANCE FOR PLANETARY MISSIONS,” 
by Daniel P. Scharf, James C. Blackmore and Behcet Acik- 
mese, filed Mar. 27, 2009; 20 

U.S. Provisional Patent Application No. 61/183,508, 
entitled “A TABLE LOOKUP SOLUTION ALGORITHM 
OF THE OPTIMAL POWERED DESCENT GUIDANCE 
PROBLEM FOR PLANETARY LANDING,” by James C. 
Blackmore, and Behcet Acikmese, filed Jun. 2, 2009; and 25 

U.S. Provisional Patent Application No. 61/286,095. 
entitled “MINIMUM LANDING ERROR POWERED 
DESCENT GUIDANCE FOR MARS LANDING USING 
CONVEX OPTIMIZATION,” by James C. Blackmore, Beh- 
cet Acikmese, and Daniel P. Scharf, filed Dec. 14, 2009. 30 

STATEMENT OF RIGHTS OWNED 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 35 
provisions of Public Law 96-517 (35 U.S.C. §202) in which 
the Contractor has elected to retain title. 

BACKGROUND OF THE INVENTION 

40 

1. Field of the Invention 

The present invention relates to systems and methods for 
guiding spacecraft, and in particular to a system and method 
for guiding a spacecraft to a soft landing on a planetary 
surface. 45 

2. Description of the Related Art 

To increase the science return of future missions to Mars, 
and to enable sample return missions, the accuracy with 
which a lander can be delivered to the Martian surface must be 
improved by orders of magnitude over the current capabili- 50 
ties. Towards this goal, our prior work developed a convex 
optimization based minimum-fuel powered descent guidance 
algorithm. Here, we extend this approach to handle the case 
when no feasible trajectory to the target exists. In this case, 
our objective is to generate the minimum landing error tra- 55 
jectory, which is the trajectory that minimizes the distance to 
the prescribed target while utilizing the available fuel opti- 
mally. This problem is inherently a non-convex optimal con- 
trol problem due to a nonzero lower bound on the magnitude 
of the feasible thrust vector. We first prove that an optimal 60 
solution of a convex relaxation of the optimal control problem 
is also optimal for the original non-convex problem, which 
we refer to as the lossless convexification of the original 
non-convex problem. Then we show that the minimum land- 
ing error trajectory generation problem can be posed as a 65 
convex optimization problem, in particular as a Second-Order 
Cone Program, and solved to global optimality with deter- 


2 

mini stic convergence. This makes the approach amenable to 
onboard implementation for real-time applications. 

SUMMARY OF THE INVENTION 

To address the requirements described above, the present 
invention discloses a method and apparatus for computing a 
thrust profile to land a spacecraft at or near a surface target, 
wherein the spacecraft is subject to a non-convex thrust con- 
straint. In one embodiment, the method comprises the steps of 
computing a minimum landing error descent solution for the 
spacecraft subject to a convexified non-convex constraint, 
wherein the descent solution includes a landing error, and 
determining the thrust profile for the spacecraft subject to the 
non-convex constraint from the computed minimum landing 
error descent solution for the spacecraft subject to the con- 
vexified non-convex constraint, and commanding the space- 
craft thrusters according to the thrust profile. In a further 
embodiment, the step of determining the thrust profile for the 
spacecraft subject to the non-convex constraint from the com- 
puted minimum landing error descent solution for the space- 
craft subject to the convexified non-convex constraint com- 
prises the steps of computing a minimum fuel descent 
solution for the spacecraft subject to the convexified non- 
convex constraint and the landing error, and determining the 
thrust profile for the spacecraft subject to the non-convex 
constraint from the computed minimum fuel descent solution 
for the spacecraft subject to the convexified non-convex con- 
straint. 

In a still further embodiment, the apparatus comprises a 
processor for computing a minimum landing error descent 
solution, comprising a memory for storing instructions 
including instructions for computing a minimum landing 
error descent solution for the spacecraft subject to a convexi- 
fied non-convex constraint, wherein the descent solution 
includes a landing error, and for determining the thrust profile 
for the spacecraft subject to the non-convex constraint from 
the computed minimum landing error descent solution for the 
spacecraft subject to the convexified non-convex constraint, 
and for commanding the spacecraft thrusters according to the 
thrust profile. The apparatus further comprises thrusters for 
applying the commanded thrust to the spacecraft. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Referring now to the drawings in which like reference 
numbers represent corresponding parts throughout: 

FIG. lisa diagram illustrating glideslope constraints for a 
minimum landing errorpowered descent. The glideslope con- 
straint requires the spacecraft to remain in a cone defined by 
a minimum slope angle a. In a minimum landing error case, 
the apex of the cone coincides with the landed position of the 
spacecraft, rather than the original target. 

FIG. 2 is a three dimensional plot of an optimal trajectory 
generated using a minimum landing error targeting approach 
for Case 1 , with r 0 = [ 1 500 0 2000,] r , r 0 =[75 0 1 00] r (units in 
meters and meters/second, respectively). In this case, a fea- 
sible solution exists, so the prioritized minimum landing error 
algorithm returns the minimum fuel solution to the target. 

FIGS. 3A-3H are diagrams presenting the result of the 
minimum landing error targeting approach for the initial con- 
ditions described with respect to FIG. 2. 

FIG. 4 is a three dimensional plot of an optimal trajectory 
using a minimum landing error targeting approach for Case 2 
with r o =[1500 0 2000, ] r , r 0 =[-75 40 100] r (units in meters 
and meters/second, respectively). In this case, no feasible 
solution to the target exists, so the minimum landing error 
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solution approach returns the solution that minimizes the final 
distance to the target, which is 268 meters. 

FIGS. 5A-5H are diagrams presenting the result of the 
minimum landing error targeting approach for the initial con- 
ditions described with respect to FIG. 4. 5 

FIG. 6 is a diagram illustrating the dependence of mini- 
mum distance from the target on time of flight for a typical set 
of initial conditions. The relationship is unimodal and has an 
optimum at 56.35 seconds. The Golden Search determines 
the optimum to be 55.83 seconds, which is an error of only to 

0. 9%. 

FIG. 7 is a diagram presenting illustrative method steps that 
can be used to command the thrusters to land the spacecraft. 

FIG. 8 is a diagram presenting exemplary method steps that 
can be used to determine the thrust profile for the spacecraft 15 
subject to the non-convex constraint from the computed mini- 
mum landing error descent solution for the spacecraft subject 
to the convexified non-convex constraint. 

FIG. 9 is a diagram of an exemplary spacecraft. 

FIG. 10 is a diagram depicting the functional architecture 20 
of a representative spacecraft guidance and control system. 

FIG. 11 is a diagram depicting an exemplary computer 
system that could be used to implement elements of the 
present invention at a ground station. 

25 

DETAILED DESCRIPTION OF PREFERRED 
EMBODIMENTS 

In the following description, reference is made to the 
accompanying drawings which form a part hereof, and which 30 
is shown, by way of illustration, several embodiments of the 
present invention. It is understood that other embodiments 
may be utilized and structural changes may be made without 
departing from the scope of the present invention. 

1. Introduction 35 

(Note: This application references a number of different 

publications as indicated throughout the specification by ref- 
erence numbers enclosed in brackets, e.g., [x] or the words 
“Ref.” or “Refs.”. A list of these different publications ordered 
according to these reference numbers can be found below in 40 
the section entitled “References.” Each of these publications 
is incorporated by reference herein.) 

The science return of previous missions to the surface of 
Mars has been limited by the accuracy with which a lander 
can be delivered to the surface. Landing accuracy is charac- 45 
terized by the 3-sigma landing ellipse, which defines the 
region around the target in which landing can be guaranteed. 
The size of this ellipse (major axis) was approximately 150 
km for Mars Pathfinder and 35 km for the Mars Exploration 
Rovers. The 2009 Mars Science Laboratory mission aims to 50 
achieve a landing ellipse of around 10 km. This means that 
landing site selection is driven by the need to find a safe 
landing site, rather than by science goals. In other words, 
landing sites must be large, flat, and relatively rock- and 
crater-free areas to ensure safe landing. These regions are 55 
usually not the sites with the maximum science return. For 
example, in the recent Phoenix mission, orbital images taken 
months before the landing showed a higher-than-expected 
concentration of large rocks at the primary landing site. This 
required a switch to an alternate landing site, which may not 60 
have been necessary if a more accurate landing were possible. 

Accurate landing is also an enabling technology for Mars 
sample return missions. One concept for such a mission 
would involve a first spacecraft gathering samples, with a 
second spacecraft arriving at Mars later to return the sample 65 
to Earth. For this to be possible, the second spacecraft has to 
land close to the first in order to enable transfer of the sample 
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without prohibitively high power and mobility requirements. 
Achieving this with a sufficiently low level of risk to the 
mission will require dramatic improvements in landing accu- 
racy. 

Recent work has focused on achieving pinpoint landing, 
which is defined as landing to within hundreds of meters of a 
target. The pinpoint landing concept consists of an entry 
phase through the Martian atmosphere, during which bank 
angle control is used to taiget the spacecraft, followed by 
parachute deployment. During the parachute phase, errors 
accumulate due to winds and atmospheric uncertainty. Once 
the parachute is released, the final powered descent phase 
then uses thrusters to land safely at the target. To do so, the 
lander needs to calculate onboard a trajectory from its a priori 
unknown location at parachute cutoff to the target. This pow- 
ered descent guidance (PDG) problem is challenging for a 
number of reasons. First, we must guarantee that any feasible 
solution obeys hard constraints, including minimum and 
maximum thrust magnitudes and a minimum glideslope 
angle. The latter constraint also prevents subsurface flight. 
Second, we must guarantee that a feasible solution will be 
found in a matter of seconds. This requirement is derived 
from the duration of the parachute and PDG phases; if the 
PDG algorithm takes too long to find a feasible solution, the 
lander can crash into the surface. Third, to both maximize the 
distance from which the lander can reach the target and mini- 
mize the amount of fuel that must be carried onboard, the 
algorithm should find the globally optimal solution. 

A great deal of prior work has developed approximate 
solutions to the powered descent guidance problem Ref. 2, 
and Refs. 5-10. The convex optimization approach poses the 
problem of minimum-fuel powered descent guidance as a 
Second-Order Cone Program (SOCP). This optimization 
problem can be solved in polynomial time using existing 
algorithms that have a deterministic stopping criterion given 
a prescribed level of accuracy. That is, the global optimum 
can be computed to any given accuracy with a deterministic 
upper bound on the number of iterations required for conver- 
gence. This is in contrast with other approaches that either 
compute a closed-form solution by ignoring the constraints of 
the problem, propose solving a nonlinear optimization 
onboard, or solve a related problem that does not minimize 
fuel use. The closed-form solution approach results in solu- 
tions that do not obey the constraints inherent in the problem, 
such as no subsurface flight constraints; in practice this 
reduces the size of the region from which return to the taiget 
is possible by a factor of ten or more. 

Nonlinear optimization approaches, on the other hand, 
cannot provide detemiinistic guarantees on how many itera- 
tions will be required to find a feasible trajectory, and are not 
guaranteed to find the global optimum. This limits their rel- 
evance to onboard applications. For more extensive compari- 
sons of the convex optimization approach to alternative 
approaches, see Refs. 10 and 14. 

We extend the convex optimization approach of Ref. 2 to 
handle the case when no feasible traj ectory to the target exists . 
If no feasible trajectory to the target exists, the onboard guid- 
ance algorithm must ensure that a safe landing occurs as close 
as possible to the original target, if necessary using all of the 
fuel available for powered descent guidance. In this paper we 
present an algorithm that solves this minimum landing error 
problem. The algorithm calculates the minimum-fuel trajec- 
tory to the target if one exists, and calculates the trajectory 
that minimizes the landing error if no solution to the taiget 
exists. In the spirit of Ref. 2, the approach poses the problem 
as two Second-Order Cone Programs, which can be solved to 
global optimality with detemiinistic bounds on the number of 
iterations required. This minimum landing error capability 
will be necessary for missions that want to increase landing 
accuracy, but cannot carry enough fuel to ensure that the 
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target can be reached in all possible scenarios. In addition the 
capability will enable safe landing as a contingency in the 
case where when the position and velocity of the lander at 
parachute cutoff is outside of the range of values anticipated 
during mission design. Causes of this misprediction include 5 
higher-than-expected winds or atmospheric parameters, such 
as atmospheric density, being outside of the modeled values. 
Since there remains a great deal of uncertainty in atmospheric 
models, this possibility should be considered in the mission 
design. to 

The key difficulty in posing the minimum landing error 
problem as a convex optimization problem is the presence of 
non-convex constraints, namely the nonlinear system dynam- 
ics and the non-convex thrust constraints. The thrust con- 
straints are non-convex because of a minimum thrust magni- 15 
tude constraint, which arises because, once started, the 
thrusters cannot be throttled below a certain level. In Ref. 2, 
we rendered the minimum-fuel guidance problem convex by 
relaxing the non-convex constraints and then proving that any 
optimal solution of the relaxed problem is also an optimal 20 
solution to the full non-convex problem. 

We refer to this approach as lossless convexification. The 
primary theoretic contribution of the present paper is an 
analogous convexification for the minimum landing error 
powered descent guidance problem. We pose two optimiza- 25 
tion problems in which the non-convex constraints have been 
relaxed to yield Second-Order Cone Programs. We solve 
these problems sequentially, and then provide a new analytic 
result that shows that any optimal solution to the second 
relaxed problem is also an optimal solution to the non-convex 30 
minimum landing error problem. This result is then used to 
develop the algorithm that solves the minimum landing error 
problem. Furthermore, we show that the solution uses the 
least fuel of all the optimal solutions to the minimum landing 
error problem. 

This paper is organized as follows: In Section II we review 35 
the convex optimization approach of Ref. 2 to solving the 
minimum-fuel PDG problem before defining the minimum 
landing error PDG problem in Section III. In Section B we 
present the main analytic result of the paper, namely the 
convexification of non-convex control constraints. In Section 40 
IV we present the new algorithm for minimum landing error 
PDG. Section V gives simulation results and Section VI pre- 
sents our conclusions. 

II. Minimum-Fuel Powered Descent Guidance 

In this section we review the minimum-fuel powered 45 
descent guidance approach of Ref. 2. The minimum-fuel 
PDG problem consists of finding the thrust profile that takes 
the lander from an initial position and an initial velocity to rest 
at the prescribed target location, and minimizes the fuel mass 
consumed in doing so. 50 

The minimum-fuel PDG problem is defined formally in 
Problem 1. Throughout this paper we use e ; . to denote a col- 
umn vector of all zeros except the \ ,h row, which is unity. We 
use to denote the two-norm of the vector. We use a.e. to mean 
almost everywhere, that is, everywhere except on a set of 
measure zero. 


(Non-convex minimum fuel guidance problem). 

max m(tf) — m(0) = min ( ^ a\\T c (t)\\dt 
tfJd) tf,T c ()J 0 

subject to: 

r(r) = g + T c (t)/m(t) m(t) = -o||T c (t)|| 

0 <Pi < IITWII sp : 


Problem 1 

(1) 60 

(2) 


-continued 


lit) etVre [0, /, ! 

(4) 

m(0) = m„„, m(tf) > m dry 

(5) 

H 0) = r„, m = h> 

(6) 

ejr(tf) = 0, [e 2 e 3 ] r r(t f ) = q, r(t f ) = 0 

(7) 


Here, q£ER 2 defines the location of the landing target on the 
surface. We use X to define the set of feasible positions of the 
spacecraft, namely the glide slope constraint: 

X={*E K R 3 :\\S(r-r(tj))\\-c T (r-r(t f ))S0}, (8a) 

where S and c are defined by the user to specify a feasible cone 
with a vertex at 


0 l 
0 0 1 


c = eytana a > 0. 


Here a is the minimum glide slope angle, illustrated in FIG. 
1 . The glide slope constraint (8a) ensures that the trajectory to 
the target cannot be too shallow and cannot go subsurface. X 
is a convex set with an interior defined by 

int A2 ={.vfE.V:3r >0 such that ctEXif\\x-c\\<r}, (9) 

and the boundary of X is given by 

5X:={xEX:x$.m\.X\. 

The inequalities in (8) ensure that the planned trajectory 
does not fly below the surface and that general cone con- 
straints on the lander’s position are satisfied. These cone 
constraints are used to specify, for example, minimum glides- 
lope angles; these ensure that the trajectory to the target 
cannot be too shallow. Equation (5) defines the initial mass of 
the lander and ensures that no more fuel than available is used. 
Equation (6) defines the initial position and velocity of the 
lander, while (7) constrains the lander to be at rest at the target 
at the final time. 

A key challenge in solving Problem 1 consists of the non- 
convex thrust magnitude constraints in (3). These prevent the 
direct use of convex optimization techniques in solving this 
problem. The key theoretical innovation of Ref. 2 is to relax 
these non-convex constraints to give the following problem 
and to show that the optimal solution of this relaxed problem 
is also an optimal solution of Problem 1 : 


(Relaxed minimum fuel guidance problem). 

Problem 2 

min | ' 1 \t) dl 

roJ„ 

(10) 

subject to: 


r(f) = g + T c (t)lm(t) m(t ) = -a-r(r) 

(11) 

IITWII Hr) 0 <pi < T(r) < p 2 

(12) 

r(r) e XV re [0, i : ] 

(13) 

m( 0) = m w „, m(tf) > m&y 

(14) 

r(0) = r 0 , HO) = r 0 

(15) 

e\r(t[) = 0, [e 2 e 3 ] r r(t f ) = q, r(tf) = 0 

(16) 


Note that the non-convex thrust constraints in Equation (3 ) 
have been replaced with convex set of constraints (Eq. 12). In 
Ref. 2 we showed that this constraint relaxation allows the 
discrete-time form of Problem 2 to be posed as a convex 
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optimization problem. Furthermore, the following lemma 
formally states that an optimal solution to the relaxed mini- 
mum-fuel problem is also an optimal solution to the full, 
non-convex minimum-fuel PDG problem. 

Lemma 1. 5 

Let {t/,T c (*),r*(t)} be an optimal solution to Problem 2 
such that the corresponding state trajectory r*(t) is on the 
boundary of the state constraints 3X for at most one point on 
te[0,t/]. Then, {ty l ',T c (*)J is an optimal solution to Problem 
1 . 10 

Proof of Lemma 1 is provided in the documents: Acikmese, 

B. and Ploen, S. R., “Convex Programming Approach to 
Powered Descent Guidance for Mars Landing,” AIAA Jour- 
nal of Guidance, Control and Dynamics, Vol. 30, No. 5; 2007, 15 
pp. 1353-1366 and Blackmore, L., and Acikmese, B., “Mini- 
mum Landing Error Powered Descent Guidance for Mars 
Landing using Convex Optimization” Under Review. Journal 
of Guidance, Control, and Dynamics, 2009, which are hereby 
incorporated by reference herein. Specifically, “Convex Pro- 20 
gramming Approach to Powered Descent Guidance for Mars 
Landing,” provides a proof of Lemma 1 when the optimal 
state trajectory is strictly in the interior of the state con- 
straints, while “Minimum Landing Error Powered Descent 
Guidance for Mars Landing using Convex Optimization” 25 
extends this to cases where there is contact with the boundary 
of the state constraints. 

The approach of Ref. 2 solves the non-convex minimum- 
fuel PDG problem by solving a relaxed, convex version of the 
problem, which is a convex optimization problem. The con- 11 1 
vexification is lossless, in the sense that no part of the feasible 
space of the original problem is removed in the convexifica- 
tion process. The resulting optimization problem is a Second- 
Order Cone Program, for which solution techniques exist that 
guarantee the globally optimal solution can be found to a 
certain accuracy within a deterministic number of iterations. 
However, if a feasible solution does not exist, the optimiza- 
tion will simply report “infeasible,” even though it may still 
be possible to land safely at some distance from the original 
target. In the next sections, therefore, we extend the approach 
of Ref. 2 to solve the minimum landing error problem. 

III. Minimum Landing Error Powered Descent Guidance 
Problem Statement 

FIG. lisa diagram illustrating glideslope constraints for a 
minimum landing error powered descent. The glideslope con- 
straint requires the spacecraft to remain in a cone 102 defined 
by a minimum slope angle a. In a minimum landing error 
case, the apex 104 of the cone 102 coincides with the landed 
position 104 of the spacecraft, rather than the original target 
106 . 

The minimum landing error powered descent guidance 
problem consists of minimizing the final distance from the 
target subject to non-convex thrust magnitude constraints and 
glideslope constraints, while ensuring that no more fuel mass 
is used than is available. The problem is stated formally in 
Problem 3. below 


(Non-convex minimum 
landing error guidance problem). 

Problem 3 

min ||r(r / )|| 2 

(17) 



subject to: 


fit) =g + T c (t)/m{t) m(t) = -ar||r c (r)|| 

(18) 
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-continued 


O 

A 

£ 

IA 

IA 

$ 

(19) 

r(r) e XV re [0, r f ] 

(20) 

m( 0) = m„„, m[tf) > m dry 

(21) 

r( 0) = r-o, KO) = f 0 

(22) 

o 

II 

o' 

II 

3 

(23) 


There are a number of key differences between this and the 
minimum-fuel guidance problem (Problem 1). First, the cost 
in Problem 3 is now the squared Euclidean distance from the 
target at the final time. Minimizing the squared distance is 
equivalent to minimizing the distance itself, since IKV)II=° 
andx 2 is monotonic forxSO.}. To simplify notation we have 
assumed that the target is at zero, without loss of generality. 
Second, the final position is no longer required to be at the 
goal as in (7). Instead, (23) constrains the final altitude to be 
zero and the final velocity to be zero. Note that the cone 
constraints in (8) are defined around the state at the final time 
step, and not around the origin. This allows glideslope con- 
straints to be imposed even in the case that it is not possible to 
reach the target, as illustrated in FIG. 1. 

Once again, the key difficulty in solving Problem 3 is that 
the thrust constraints (19) are non-convex. The approach of 
Ref. 2 suggests that we overcome this problem through a 
similar convexification of the thrust constraints. Unfortu- 
nately, the result in Lemma 1 does not apply to Problem 3. 
This means that we are not guaranteed that an optimal solu- 
tion to a convexified relaxation of the minimum landing error 
problem is an optimal solution to the minimum landing error 
problem. In Section IV we therefore propose a new convexi- 
fication for minimum landing error powered descent guid- 
ance. 

IV. Minimum Landing Error Powered Descent Guidance 
Technique 

In this section we present a new technique for minimum 
landing error powered descent guidance. In Section W-A we 
present the algorithm. In Section IV-B we present a corre- 
sponding proof that the algorithm returns the globally optimal 
solution to Problem 3 if a solution exists. 

A. Description 

In this section we describe the main algorithm for the 
solution of the minimum landing error powered descent prob- 
lem (Problem 3). The key idea is to perform a prioritized 
optimization such that, first, the distance between the pre- 
scribed target and the achievable landing location is mini- 
mized by solving Problem 4. Then a minimum fuel trajectory 
achieving this minimum landing error is generated by solving 
a slightly modified version of Problem 5. As we show in 
Section IV-B, this approach ensures that the resulting trajec- 
tory satisfies non-convex thrust constraints and returns the 
globally optimal solution to Problem 3 if one exists. The new 
algorithm is given in Table 1 . 

TABLE 1 


Prioritized Powered Descent Guidance Algorithm 
Prioritized Powered Descent Guidance Algorithm 


1) Solve the relaxed minimum landing error guidance problem (Problem 
4) for {ty*, T c *(-), r*(-)} with corresponding trajectory r*(-). If no 
solution exists, return infeasi le 

2) Solve the relaxed minimum- fuel guidance problem to specified range 
(Problem 5) with D = ||r(t/)|| for {t/, T/Q, T 1 ^)}. 

3) Return {t/, T/(-)}. 
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(Relaxed minimum landing error guidance problem). 

Problem 4 

min ||r(r / )|| 2 
t f ,T c (), r() 

(24) 

subject to: 


r(r) = g + T c (t)/m(t ) fn{t) = — orr(r) 

(25) 

IITcWII S r(f) 0 < pi < r (r) <p 2 

(26) 

rU ) e XV te [0, ly ] 

(27) 

m(0) = m w „, mitf) > m drf 

(28) 

r( 0 ) = r 0 , r(0) = r 0 

(29) 

r(tf) T e i = 0, r(tf) = 0 

(30) 


(Relaxed minimum fuel 

guidance problem to specified range). 

Problem 5 

min i ^ r(t)dr 

< f .TM ro J 0 

(31) 

subject to: 


r(t) = g + T c (t)/m(t) m(t) = -aT(r) 

(32) 

IILWII £ T(r) 0<p! <T(r)<p2 

(33) 

r(t) r e i > 0 V r e [0, tj] 

(34) 

r(r) e XV t e [0, t { ] 

(35) 

m(0) = m„„, m(t f ) > m dr) 

(36) 

r(0) = r 0 , r(0) = r 0 

(37) 

r(t f ) r e i = 0, ||r(r/)|| < D, r(t f ) = 0 

(38) 


To generate Problem 4 we relaxed the non-convex thrust 
constraints of the original minimum landing-error problem 
stated in Problem 3. This relaxation is performed in the same 
manner that the minimum-fuel PDG problem is relaxed from 
Problem 1 to Problem 2. Since the relaxed minimum landing- 
error problem (Problem 4) has convex inequality constraints 
on the state as well as on the controls, it can be solved with 
existing solvers. However, unlike minimum-fuel PDG, in 
which an optimal solution to the relaxed Problem 2 is also an 
optimal, feasible solution to Problem 1 , an optimal solution to 
the relaxed minimum landing-error problem (Problem 4) is 
not necessarily an optimal, feasible solution to the original 
minimum landing-error problem, Problem 3. In particular, 
another step is needed to ensure that the non-convex thrust 
constraints are satisfied. This step consists of solving Problem 
5, which also has relaxed constraints. However, since Prob- 
lem 5 minimizes fuel use, we can use the result of Lemma 1 
to prove that the solution to Problem 5 will satisfy the non- 
convex thrust constraints of Problem 3. We will show in 
Section IV-B that the solution to Problem 5 is the optimal 
solution to Problem 3, and that it uses the minimum possible 
fuel of all optimal solutions to Problem 3. 

In order to solve Problems 4 and 5 in practice, three addi- 
tional steps are required; first a change of variables to remove 
the nonlinear (and hence non-convex) dynamic constraints; 
second, a discretization in time; and third, a line search for the 
optimal time of flight. For these steps we use an identical 
approach to Ref. 2, which we review in Sections IV-C through 
IV-E. 


10 

B. Analytic Convexification 

In this section we provide two necessary technical lemmas, 
before proving the main convexification result, which is given 
in Theorem 1 . 

5 Le mm a 2: 

Assume that an optimal solution to Problem 5 exists, which 
we denote {t/ ,T/(*),r 5 (*)}, which has a corresponding tra- 
jectory r 5 ^) that is on the boundary 3X for at most one point 
on t<E[0, t^|. Then {t/,T/(*)} is a feasible solution to prob- 
10 lem 1 with q=[e 2 e 3 ] r(t/). 

Proof: 

We first claim that {t / - 5 ,T c 5 ( l P),r 5 ( l F)} is an optimal solu- 
tion to problem 2 with q=[e 2 e 3 ] r r 5 (t/). The proof of this is by 
contradiction. Let us assume that there exists a solution {y*, 
15 T c *(*),r*(*)} that satisfies the constraints of Problem 2 with 
q=[e 2 e 3 ] 7 r 5 (L 5 ), but that also lias: 


C'f r „ . r/ r5 , . (39) 

r (r)dt< r \t)dr. 

2U JO JO 

Comparing constraints, we see that {y*,T t ,*( , ),r*(»)} is a 
feasible solution to Problem 5. Hence {y*,T £ .*( , ),r*(*)} is 
, . also a feasible solution with lower cost than the optimal 
solution {t/,T/(*),r 5 (*)}, which leads to a contradiction. 
Hence {W(*)U( •)} is an optimal solution to Problem 2 
with q=[e 2 e^r^y 5 ). Since r 5 (t) is on the boundary 3X for at 
most one point on t£[0,y 5 ], then by Lemma 1, {y 5 ,T c 5 (*),T 5 
30 (*)} is a feasible solution to Problem 1 with q=[e 2 e 3 ] r r 5 (t / - 5 ). 

Lemma 3: 

Assume that an optimal solution to Problem 4 exists, which 
we denote {t/,Tj(.),U( •)}, which has a corresponding tra- 
jectory r 4 («). Then there exists a feasible, optimal solution to 
35 Problem 5 with D=||r 4 (t / 4 )||, which we denote {t/,T/(*),r 5 
(•)} with corresponding trajectory r 5 ^). Assume further that 
this trajectory is on the boundary 3X for at most one point on 
t£[0,y 5 ]. Then {t/,T c 3 (*)} is a feasible, optimal solution to 
Problem 3. Furthermore {y 5 ,T e 5 (*)} uses the minimum pos- 
40 sible fiiel of all optimal solutions to Problem 3. 

Proof: 

First, note that {t / ?*,T t . 4 *(*),r 4 (')} is a feasible solution to 
Problem 5 with D=||r 4 (y 4 )|| since the only additional con- 
straint in Problem 5 is l|r(y 4 )ll SD, and by design, we have 
D=||r 4 (y 4 )||. Hence, a feasible solution to problem 5 exists for 
D=jjr 4 (t / - 4 )jj. Solving this problem, we obtain {y 5 ,T c 5 (*),T 5 
(•)}, and from Lemma 2, we know that {y 5 ,T c 5 (*)} is a fea- 
sible solution to Problem 1 for q=[e 2 e 3 ] r r(y 5 ). By compar- 
ing constraints, we see that any feasible solution to Problem 1 
5Q is a feasible solution to Problem 3. Hence, {y 5 ,T c 5 (*)} is a 
feasible solution to Problem 3. This means that: 
||^(r/)||S||H(r/)||, 

Where r 3 !*) denotes the trajectory corresponding to the 
optimal solution to Problem 3. Now since Problem 4 is a 
55 relaxation of Problem 3, we know that ||r 4 (y 4 )||£=||r 3 (y 3 )||. 
Since in Problem 5 we have assigned D=||r 4 (t / 7)||, we also 
know that ||r 5 (y 5 )||S||r 4 (y 4 )|| and hence: 

||r 5 (f/)[K||^(f/)||, 

60 Combining (40) and (41) we have ||r 5 (t / 5 )||S||r 3 (t / 3 )||. Hence 
the landing error in the optimal solution to Problem 5 is the 
same as that m a globally optimal solution to Problem 3. We 
have already shown that {y?,T/(*)} is a feasible solution to 
Problem 3. Hence {ty 5 ,'l) 5 (*)J is an optimal solution to Prob- 
65 lem 3. 

We now show that {y 5 ,T/(*)} uses the minimum possible 
fuel of all optimal solutions to Problem 3. We know that 
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{ty 5 ,T £ , 5 (*)} is ail optimal solution to Problem 3. Hence all 
optimal solutions to Problem 3 have |[r(t / )||=||r 5 (t / . 5 )||S||r 4 (t / 4 )||. 
The proof is by contradiction. Assume that there exists an 
optimal solution to Problem 3 with corresponding 

trajectory r + (*) that also has: 


t 5 

pj|7?WII<lr< pj|r c 5 WII dr- 
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cal sense. Our approach removes this ambiguity from the 
problem description and obtains a physically meaningful 
solution. 

C. Change of Variables 

In this section we review the change of variables employed 
by Ref. 2 to remove the non-convex constraints introduced by 
the nonlinear state dynamics (25). The change of variables 
leads to Problem 6, presented later in this section, that is a 
continuous-time optimal control problem with a convex cost 
and convex constraints. The change of variables is given by: 


Comparing constraints, since ||r f (t/)||=||r 4 (U)|| we know that 
{t / t ,T c , 1 '(*),r 1 '(*)} with r t (*)=||T £ t (*)|| is also a feasible solu- 
tion to Problem 5 with D=||r 4 (t/||. From (42) this solution has 
a lower cost than {t / - 5 ,T ( . 5 (*)}. This leads to a contradiction 15 
since {t/,T c 5 (*),T5 3 (*)} is the optimal solution to Problem 5 
with D=||r 4 (t/||. Hence there is no optimal solution to Problem 
3 that uses less fuel than {t / r 5 ,T/(*)}. 

The following theorem is the main result of this paper and 
it follows from Lemmas 1 through 3. 20 

Theorem 1 . If a solution to the non-convex minimum land- 
ing error guidance problem (Problem 3) exists, then the pri- 
oritized powered descent guidance algorithm returns a solu- 
tion {t^jT^*)} with trajectory r t (*). If this trajectory is on the 
boundary of the state constraints 0X for at most one point on 25 
the interval [0,ty ], then it is an optimal solution to Problem 3. 
Furthermore the returned solution uses the minimal fuel 
among all optimal solutions of Problem 3. 

Proof: Since Problem 4 is a relaxation of Problem 3, we 
know that if there is a feasible solution to Problem 3, there 30 
exists a feasible solution to Problem 4, and the prioritized 
powered descent guidance will not return infeasible. Then 
from Lemma 3 we know that {ty ,T,2(*)} is an optimal solu- 
tion to Problem 3, and that this solution uses the least fuel of 
all optimal solutions to Problem 3. 35 

Note that Theorem 1 implies that the prioritized powered 
descent guidance algorithm returns infeasible only if no solu- 
tion to the non-convex minimum landing error guidance 
problem (Problem 3) exists. Hence the convexification 
approach is lossless, in the sense that no part of the feasible 40 
region of the original problem is removed by convexifying the 
non-convex constraints. 

Remark 1. As with the minimum fuel powered descent 
guidance problem we have observed that, for Mars landing, 
all the optimal trajectories that are generated via solving the 45 
relaxed minimum-fuel guidance problem to specified range 
touch the boundary of the feasible state region at most one 
time instant. 

Hence the prioritized powered descent guidance algorithm 
has always returned optimal solutions to the original non- 50 
convex minimum landing error problem for Mars powered 
descent guidance. This includes an extensive empirical inves- 
tigation across the space of feasible initial conditions and 
system parameters. 

Remark 2. Observe that in Step 1 ofthe prioritized powered 55 
descent guidance algorithm, we do not necessarily obtain a 
fuel-optimal solution, but that Step 2 ensures a fuel-optimal 
solution is found. In this way the two-step approach is a way 
of prioritizing a multi-objective optimization problem. This 
approach is different from the more typical regularization 60 
procedure where both distance and fuel costs are combined in 
a single cost function to ensure that a single optimal solution 
exists. The prioritization approach removes two of the key 
problems associated with regularization. First, regularization 
requires that the relative weights on fuel and distance are 65 
chosen, which is usually carried out in an ad-hoc manner. 
Second, regularized solutions do not necessarily make physi- 


a r 

<T = — 

111 

A T c 
u — — 
m 

A i 

I = In m. 

Equation (25) can then be rewritten as: 


(43) 


f(r) = u(t) + g. 


(44) 


m(r) 

i = —— = - acr(t ). 
m(t) 


(45) 


The change of variables therefore yields a set of linear 
equations for the state dynamics. The control constraints, 
however are no longer convex. These are now given by: 


v*e[o ,t f ] 

(46) 

p i e~ z(l) ^o(t) ^ p 2 e ~ z< ' l> V?E/0 ,tj\. 

(47) 


The approach of Ref. 2 uses a second order cone approxi- 
mation of the inequalities in (47) that can be readily incorpo- 
rated into the SOCP solution framework. The left-hand 
inequality of (47) can be approximated by a second order 
cone by using the first three terms of the Taylor series expan- 
sion of e _z , giving: 


, (z-zo) 2 l 

pie 41 p- (z-zo) + - J <o% 

where z 0 is a given constant. For the right-hand inequality in 
(47), a linear approximation of e _z is used corresponding to 
the first two terms of the Taylor series expansion of e~ z , 
thereby obtaining: 

agpjfT^l-fe-Zo)]. 

A linear approximation must be used, since requiring a 
variable to be less than a quadratic is non-convex. Letting: 


pi — pi e 'll, p 2 =p 2 c 


(48) 


we obtain the following second order cone and linear approxi- 
mations of (47): 


piW 1 - (z(t) - zo(t)) + 


V r e [ 0, tf]. 


fe(0 - zoWr 


(49) 


ait) < p 2 (r)[l - (z(t) - Zo(0)l 
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where: 


z 0 (7)=ln(m we/ -ap 2 ?), (50) 

and m, a is the initial mass of the spacecraft. An approxima- 
tion of Problem 4 can now be expressed in terms of the new 
control variables: 



t e [0, t f ]. 


(60) 


(Relaxed minimum landing error 
problem with changed variables). 

Problem 6 

min ||r(r / )|| 2 
tf,T c {\ no J 

(51) 

subject to: 


r(t) = g + u(t) t(t) = -acr(t) 

(52) 

ll«WII T(r) 

(53) 

, , , (z(t) - zo(0 ) 2 1 

(54) 

- (z(r) - zo(t)) + ^ - 


o-(f) < (t 2 (r)[l - (z(t) - ZoW)l 


r(t) eXVts [0, iy] 

(55) 

m(0) = m„„, m(tj) > m dry 

(56) 

r(0) = r 0 , 7(0) = 7 o 

(57) 

r(t f fei = 0, r(t f ) = 0. 

(58) 


Problem 6 is an approximation of the relaxed minimum 30 
landing error PDG problem (Problem 4) in which the nonlin- 
ear equality constraints have been eliminated. Furthermore, 
Ref. 2 shows that the approximation of the inequalities in (47) 
given by (49) is generally very accurate for both parts of the 
inequality, and derives an analytic upper bound on the 35 
approximation error. Problem 6 is a continuous-time optimal 
control problem with convex constraints and a convex cost 
function. To solve this using a direct optimization approach, a 
discretization in time is required, which we describe in Sec- 
tion IV-D. 40 

D. Time Discretization 

In this section, we apply the discretization of Ref. 2 to 
Problem 6, and develop a numerical algorithm to solve the 
resulting discrete version of the problem. The discretization 
of Problem 6 converts the infinite-dimensional optimization 45 
problem to a finite-dimensional one by discretizing the time 
domain into equal time intervals and imposing the constraints 
at edges of the timesteps, which we refer to as temporal nodes. 
Since the constraints are linear or second order cone con- 
straints, the resulting problem is a finite-dimensional SOCP 50 
problem that can be efficiently solved by readily available 
algorithms. 

For any given time interval [0,t^, and time increment. At, 
the temporal nodes are given as: 


t k =kAt k= 0, . . . , N, where NAt=tf. Define a vector of 
parameters as: 


Then the solution of the differential equations (44) and (45) at 
the temporal nodes and the control inputs at the temporal 
nodes can be expressed in terms of these coefficients: 



>(!*) 1 


ro 

A 

yt = 

7(1*) 

= 4>i 

>0 


z(ft)J 


. In m wet 

k = 1 

, ... ,N 



.<rt 


A 


u(h) 
oih ) . 


= T tn 


k = l N. 


(61) 


(62) 


where A k and Y^. are matrix functions of the time index 

k determined by the basis functions chosen. In this paper, we 
use piecewise linear basis functions with M=N, such that: 


m = 


4^2 when t e [r y— i , tj) 

t-tj 

- when I e [tj, tj+i) 
0 otherwise. 


and: 


Pj = 


u(tj) 

°itj) 


j = o M. 


(63) 


(64) 


This corresponds to first-order hold discretization of a lin- 
ear time invariant system for the dynamics of the spacecraft 
(44) and (45), with the vector [r(\ k ) T ,r(\ k ) T ,z(\. k ) T ] T &$ the state. 
Explicit expressions for <S> k , 'P Ar , A k and Y k can be obtained 
using standard techniques; we do not repeat this derivation 
here. As noted by Ref. 2, more sophisticated basis functions 
such as Chebyshev polynomials can be used, which may 
allow the controls to be described with significantly fewer 
coefficients, i.e. M«N. Now, with the following additional 
notation, Problem 7 describes the discretized version of Prob- 
lem 6: 

e=P 3x3 o 3x JF=io 1x6 1 ] £„ = 7 / 3,3 0 3x JE=[0 3x3 

/jx3 63, ,] - ( 65 ) 


Po 


n = 


Pm 


(59) 60 


(Discretized Relaxed Minimum 
Landing Error Guidance Problem). 

min||£>,v || 2 

N,tj 


where p^E ® 4 . We describe the control input, u, and a in terms 65 
of these parameters and some prescribed basis functions, 

<h(‘), • • • , <M*) : 


subject to: 

\\E u Y t i 1 \\<elY kn k = 0,... ,N 


Problem 7 


(66) 


(67) 
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w(f) U 


-continued 


~(Fy k -zoUk)) + 


( Fy k -zo(tk)) 2 
2 


elYtri < )/ 2 (r)[ 1 - ( Fy t -zo(tO)l 


Ey k eXk = 1, ... , N 
Fy N > lum,^ 

yJi e i = o, = 0 


y* =4>* 


''o 

h) 

trim,. 


+ V t >7* = l, 


( 68 ) 


(69) 

(70) 

(71) 

(72) 


Note that, for any given N, Problem 7 defines a finite- 
dimensional second-order-cone program (SOCP), which can 
be solved very efficiently with guaranteed convergence to the 
globally optimal solution by using existing SOCP algorithms. 

Here N describes the time of flight since t ,=NAt. To find the 20 
optimal time of flight, Ref. 2 propose performing a line search 
for the optimal N, solving at each iteration an SOCP for the 
remaining optimization parameters T). We perform an identi- 
cal search to solve Problem 7, as described in Section IV-E. , . 

E. Time of Flight Search 

For minimum-fuel powered descent guidance, Ref. 2 uses 
a line search to find the optimal time of flight X*. In this 
section we apply this approach to the minimum landing error 
guidance problem. Extending Ref. 2, Ref. 1 9 uses a Golden 30 
Search technique, which ensures that the interval in which the 
optimal value is known to lie shrinks by the same constant 
proportion at each iteration. Golden Search has been shown to 
be robust and efficient, and it gives an explicit interval in 
which the optimum is known to lie. This last property means 35 
that the search can be terminated when sufficient accuracy has 
been achieved. The technique relies, however, on two key 
properties of J(y), the optimal cost of Problem 3 as a function 
of U 

J AC\ 


J(t f ) = 


min 

T c an) 


||/-(f / )|| 2 subject to (18) through (23). 


(73) 


approach checks analytically whether t„ is a true upper bound. 
If not, t„ is increased until it is an upper bound on X*. 

In Ref. 19, the authors observed experimentally that is 
indeed a unimodal function. Since the present paper is con- 
cerned with the minimum landing error problem, this conclu- 
sion does not carry over from the minimum-fuel case. In 
Section V-B we show experimentally that J(t) is unimodal, 
and that the Golden Search approach finds t^to within a few 
percent. 

V. Simulation Results 

In this section we present simulation results obtained using 
the new algorithm. The Second Order Cone Programs were 
solved using G-OPT, a JPL in-house convex optimizer. Simu- 
lations were run on a Macbook Pro 2.4 GHz with 4 GB RAM. 
In Section V-A we present some example solutions generated 
by the new approach, while in Section V-B we demonstrate 
empirically the unimodality of the optimal solution with 
respect to the time of flight. 

A. Example Solutions 

We first consider a case where the target \emph{can} be 
reached given the available fuel, then a case when the taiget 
\emph{cannot} be reached. The spacecraft parameters for 
these simulations are: 

g=[- 3.71 14 0 0] r Wj o ,=1.505 kg m w ,=1905 kg 

I sp =225 s pi=4972N p 2 =13260N (74) 

A glideslope constraint is used to prevent the trajectory 
from entering at a more shallow angle than 4°. The spacecraft 
initial position is given by: 


r o = 


1500 m 
0 

2000 m 


(75) 


In Case 1. the initial velocity is denoted r 0 (1) and is given 
by: 


-75 m/s 
0 

100 m/s 


(76) 


First, we must know an interval in which X* is known to lie. 
That is, we must find bounds X, and t u such that E=ty !, =t u . 
Second, J(t^) must be unimodal. 

In Ref. 19, the authors solve the one-dimensional powered 
descent guidance problem using the approach of Ref. 5 to 
obtain values for t, and t i( . In the case of minimum landing 
error guidance, we can use the same approach to obtain f. 

The approach of Ref. 5 gives the minimum-time solution in 
the vertical dimension only; that is, it determines the mini- 
mum-time thrust profile that eliminates the initial vertical 
velocity and ends with the lander at zero altitude. Denote the 
optimal time of this solution as X f ' n . In the minimum landing 
error guidance problem the lander must also have zero final 
altitude and vertical velocity, hence the constraint set for the 
minimum landing error guidance problem is tighter than for 
the one-dimensional problem. This means that any feasible 
time-of-flight for Problem 3 is no less than Xj l , and hence we 
can set Xj=X^ D ^X* . 

Hie approach of Ref. 1 9 for obtaining t„ does not, however, 
extend to the minimum landing error case. Instead, we use a 
heuristic scaling to set t u =k scaie t/, where k scale is on the order 
of 3. Hien, assuming unimodality of J(t^, the golden search 


In Ref. 2 it is shown that there is a feasible solution to the 
target in this case, but that the solution requires almost all of 
the available fuel mass. This example illustrates the value of 
the convex optimization, which guarantees finding a solution 
if one exists; even cases at the edge of the physical feasibility 
can be solved. In FIG. 2 and FIGS. 3A-3D, we show results 
generated for Case 1 using the new prioritized minimum 
landing error PDG approach (Table 1) with 55 time discreti- 
zation points. FIG. 2 is a three dimensional plot of an optimal 
trajectory 200 generated using a minimum landing error tar- 
geting approach for Case 1 , with r 0 = [ 1 500 0 2000 ,] T , r 0 =[7 5 
0 1 00 ] T (units in meters and meters/second, respectively). In 
this case, a feasible solution exists, so the prioritized mini- 
mum landing error algorithm returns the mini mum fuel solu- 
tion to the target. FIGS. 3A and 3B show the optimal trajec- 
tory 200 in horizontal and vertical plane transfers, while 
FIGS. 3C and 3D show the angle above the surface 302 and 
optimal throttle 304, respectively. FIGS. 3E-3G illustrate the 
position, velocity, control acceleration for the optimal trajec- 
tory 200 inX, Y and Z coordinates, while FIG. 3G illustrates 
the thrust angle and rate. Golden Search was terminated when 
the optimal time of flight was known to within an interval of 
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3.0 s. Since there exists in this case a feasible solution to the 
target, the algorithm returns the minimum-fuel solution to the 
target. This solution requires 399.5 kg of fuel, has ty*=78.3 s 
and is identical to the solution reported by Ref. 2. The total 
computation time required was 14.3 s, and 23 iterations of 5 
Golden Search were used. 

In Case 2, there is an additional initial velocity in the y 
direction: 


—75 m/s 
40 m/s 
100 m/s 


10 

(77) 


15 


Since Case 1 used almost all of the available fuel, and Case 
2 has an initial velocity in the y direction away from the target, 
it is most likely that there will be insufficient fuel to reach the 
target. We verily this by attempting to solve the minimum- 
fuel PDG problem using the algorithm of Ref. 2, which 
reports that the problem is infeasible, and that at least 410 kg 
of fuel is required to reach the target. The minimum landing 
error targeting algorithm, however, finds a solution that 
ensures safe landing at a distance of 268 m from the target. 
The solution, shown in FIG. 4 and FIGS. 5A-5H have has 
ty=77.6 s and uses the full 400 kg of available fuel. The total 
computation time required was 16.24 s, and 23 iterations of 
Golden Search were used. FIGS. 5 A and 5B show the result- 
ing optimal trajectory 400, while FIG. 5C and FIG. 5D show 
the angle above the surface 502 and the optimal throttle 504, 
respectively. FIGS. 5E-5G show the values of position, veloc- 
ity, and control acceleration for the trajectory as a function of 
time, respectively, while FIG. 5H shows the thrust angle and 
rate. 

B. Unimodality of Cost as a Function of Time of Flight 

In this section we demonstrate empirically that the optimal 
cost of the minimum landing error solution is unimodal in the 
time of flight. Throughout this section we use the spacecraft 
parameters given in (74). FIG. 6 shows a plot 600 J(ty) for a 
typical set of initial conditions, illustrating the dependence of 
minimum distance from the target on time of flight for a 
typical set of initial conditions. The graph was generated by 
specifying tyin increments of 1 s and solving Problem 7 for 
each value of ty In this case J(ty) is clearly unimodal, as 
required. By reducing the time increments to 0.01 s close to 
the minimum, the optimum was found to be at ^=56.35 s. The 
Golden Search determines the optimum to be 55.83 s, which 
is an error of only 0.9%. Again Golden Search was terminated 
when the optimal time of flight was known to within an 
interval of 3.0 s. The efficacy of the Golden Search approach 
was investigated for a range of initial conditions by calculat- 
ing the error between ty* determined through Golden Search, 
and the true optimum. The true optimum was determined 
approximately by calculating J^) for ^[30,150] in incre- 
ments of 1 s. Initial conditions were selected at random using 
a uniform distribution across a box of values specified by: 
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Since we are only concerned with the unimodality of J(ty) in 65 
the minimum landing error case, solutions for which a fea- 
sible trajectory to the target existed were discarded. One 


hundred random initial conditions were chosen, and for these 
the average percentage error in ty* was 0.012% with a stan- 
dard deviation of 0.057%. The maximum error across all 
instances was 1.8%. This demonstrates that the Golden 
Search approach is effective for a broad range of initial con- 
ditions in the case of minimum landing error powered decent 
guidance. 

VI. Implementation 

A. Technique 

FIG. 7 is a diagram presenting illustrative process steps that 
can be used to land the spacecraft. A minimum landing error 
descent solution is computed, as shown in block 702. As 
described in Section IV above, the solution can be computed 
for a spacecraft that is subject to non-convex thruster con- 
straints, such as a minimum thrust and a maximum thrust. 
This solution is computed, by convexifying the non-convex 
thrust constraint, and determining a solution to the system 
subject to the convexified non-convex thrust constraint, as 
described in Problem 4. 

If no computed solution exists, then landing cannot be 
performed safely under the specified conditions, and an 
“infeasible” result is returned. If a computed solution exists, 
it may be one in which the landing error range (difference 
between the desired landing point and the computed landing 
point, or D=||r(t/)||) is zero, or it may be one in which the 
landing error range D is non-zero (e.g. the minimum landing 
error solution result in a landing error). 

Next, a thrust profile for the spacecraft is determined from 
the result obtained above, as shown in block 704. In a simple 
embodiment, this can be accomplished, for example, by sim- 
ply setting the thrust profile to T c *(*) as obtained by solving 
Problem 4. This thrust profile may be used to command the 
spacecraft thrusters during powered descent, as shown in 
block 706. 

However, as described above, the optimal, feasible solution 
to the convexified minimum landing error problem (Problem 
4) is not necessarily an optimal, feasible solution to the origi- 
nal (non-convex) minimum landing error problem. Also, 
while the solution to Problem 4 computed above may provide 
a solution to land the spacecraft on the surface with minimum 
error when feasible, this is not necessarily a fuel-optimal 
solution. That is, the solution is not necessarily that solution 
that minimizes landing thruster fuel use in landing the space- 
craft at the landing point. Therefore, setting the thrust profile 
to T c *(*) for a this point will not guarantee a solution that (1) 
minimizes landing error, (2) minimizes landing thruster fuel 
use, and (3) ensures that all of the non-convex thruster con- 
straints are satisfied. 

A solution that meets all three objectives can be obtained, 
however, by taking the further step of convexifying the non- 
convex thruster constraints and solving the convexified mini- 
mum fuel guidance problem subject to constraints that 
include the landing error range D computed in block 702 
above (e.g. solving Problem 5). The resulting solution mini- 
mizes landing error and landing thruster fuel use and ensures 
that all of the non-convex thruster constraints are satisfied. 

FIG. 8 is a diagram illustrating how the process of block 
704 (determining the thrust profile for the spacecraft subject 
to the non-convex constraint from the solution computed in 
block 702) can be performed in such a way so as to minimize 
landing error and fuel use, and so as to ensure that all non- 
convex thruster constraints are satisfied. As shown in block 
802, a minimum fuel descent solution for the spacecraft is 
computed subject to the convexified non-convex constraint 
and the landing error D computed in block 702. The resulting 
thrust profile T c + (*) for the spacecraft subject to the non- 
convex thrust constraint is then computed from the minimum 
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fuel descent solution for the convexified non-convex thrust 
constraint, as shown in block 804. This resulting thrust profile 
T c + (*) is not only a minimum error solution and a minimum 
fuel solution to land the spacecraft with that minimum error, 
it is also a solution that satisfies all non-convex thruster con- 
straints are satisfied. 

B. Spacecraft and Control System 

FIG. 9 illustrates a three-axis stabilized spacecraft 900. 
The spacecraft 900 has a main body 902, a high gain narrow 
beam antennas 906, and a telemetry and command antenna 
908 which is aimed at a control ground station. The spacecraft 
900 may also include one or more sensors 912 to measure the 
attitude of the spacecraft 900. These sensors may include sun 
sensors, earth sensors, and star sensors. 

The spacecraft 900 has a power supply that may include an 
internal nuclear power supply and/or solar panels 904 that are 
used to collect solar energy. The solar panels 904 may be 
stowed or jettisoned prior to landing. 

Since the solar panels are often referred to by the designa- 
tions “North” and “South”, the solar panels in FIG. 9 are 
referred to by the numerals 904N and 904S for the “North” 
and “South” solar panels, respectively. The spacecraft may 
have a plurality of attitude control thrusters as well as landing 
thrusters 912, which are used when the spacecraft 900 is 
landing on a surface. 

The three axes of the spacecraft 900 are shown in FIG. 9. 
The pitch axis Y lies along the plane of the solar panels 904N 
and 904S. The roll axis X and yaw axis Z are perpendicular to 
the pitch axis P and lie in the directions and planes shown. 

FIG. 10 is a diagram depicting the functional architecture 
of a representative spacecraft guidance and control system. 
Control of the spacecraft is provided by a computer or space- 
craft control processor (SCP) 1002. The SCP performs a 
number of functions which may include post ejection 
sequencing, transfer orbit processing, acquisition control, 
stationkeeping control, normal mode control, mechanisms 
control, fault protection, and spacecraft systems support, 
landing control, and others. 

Input to the spacecraft control processor 1002 may come 
from a any combination of a number of spacecraft compo- 
nents and subsystems, such as a transfer orbit sun sensor 
1004, an acquisition sun sensor 1006, an inertial reference 
unit 1008, a transfer orbit sensor 1010, an operational orbit 
sensor 1012, a normal mode wide angle sun sensor 1014, a 
magnetometer 1016, and one or more star sensors 1018. 

The SCP 1002 generates control signal commands 1020 
which are directed to a command decoder unit 1022. The 
command decoder unit operates the load shedding and battery 
charging systems 1024. The command decoder unit also 
sends signals to the magnetic torque control unit (MTCU) 
1026 and the torque coil 1028. 

The SCP 1002 also sends control commands 1030 to the 
thruster valve driver unit 1032 which in turn controls the 
landing thrusters 912 and the attitude control thrusters 1036. 

Wheel torque co mm ands 1062 are generated by the SCP 
1002 and are commimicated to the wheel speed electronics 
1038 and 1040. These effect changes in the wheel speeds for 
wheels in momentum wheel assemblies 1042 and 1044, 
respectively. The speed of the wheels is also measured and fed 
back to the SCP 1002 by feedback control signal 1064. 

Hie spacecraft control processor 1002 also sends jack- 
screw drive signals 1066 to the momentum wheel assemblies 
1042 and 1044. These signals control the operation of the 
jackscrews individually and thus the amount of tilt of the 
momentum wheels. The position of the jackscrews is then fed 
back through command signal 1068 to the spacecraft control 
processor. The signals 1068 are also sent to the telemetry 
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encoder unit 1058 and in turn to the ground station 1060. The 
spacecraft control processor 1002 also sends signals to solar 
wing drives 1046 and 1048. 

The spacecraft control processor also sends command sig- 
5 nals 1054 to the telemetry encoder unit 1058 which in him 
sends feedback signals 1056 to the SCP 1002. This feedback 
loop, as with the other feedback loops to the SCP 1002 
described earlier, assist in the overall control of the space- 
craft. The SCP 1002 communicates with the telemetry 
to encoder unit 1058, which receives the signals from various 
spacecraft components and subsystems and then relays them 
to the ground station 1060. 

The wheel drive electronics 1038, 1040 receive signals 
from the SCP 1002 and control the rotational speed of the 
15 momentum wheels. The jackscrew drive signals 1066 adjust 
the orientation of the angular momentum vector of the 
momentum wheels. This accommodates varying degrees of 
attitude steering agility and accommodates movement of the 
spacecraft as required. 

20 The use of reaction wheels or equivalent internal torquers 
to control a momentum bias stabilized spacecraft allows 
inversion about yaw of the attitude at will without change to 
the attitude control. In this sense, the canting of the momen- 
tum wheel is entirely equivalent to the use of reaction wheels. 
25 The SCP 1002 may include or have access to memory 
1070, such as a random access memory (RAM). Generally, 
the SCP 1002 operates under control of an operating system 
1072 stored in the memory 1070, and interfaces with the other 
system components to accept inputs and generate outputs, 
30 including commands. Applications ranting in the SCP 1002 
access and manipulate data stored in the memory 1070. The 
spacecraft 900 may also comprise an external communication 
device such as a spacecraft link for communicating with other 
computers at, for example, a ground station. If necessary, 
35 operation instructions for new applications can be uploaded 
from ground stations. 

In one embodiment, instructions implementing the operat- 
ing system 1072, application programs, and other modules 
are tangibly embodied in a computer-readable medium, e.g., 
40 data storage device, which could include a RAM, EEPROM, 
or other memory device. Further, the operating system 1072 
and the computer program are comprised of instructions 
which, when read and executed by the SCP 1002, causes the 
spacecraft processor 1002 to perform the steps necessary to 
45 implement and/or use the present invention. Computer pro- 
gram and/or operating instructions may also be tangibly 
embodied in memory 1070 and/or data communications 
devices (e.g. other devices in the spacecraft 900 or on the 
ground), thereby making a computer program product or 
50 article of manufacture according to the invention. As such, the 
tenns “program storage device,” “article of manufacture” and 
“computer program product” as used herein are intended to 
encompass a computer program accessible from any com- 
puter readable device or media. 

55 C. Ground Station Processing 

FIG. 11 is a diagram illustrating an exemplary computer 
system 1100 that could be used to implement elements of the 
present invention at the ground station. The computer 1102 
comprises a general purpose hardware processor 11 04 A and/ 
60 or a special purpose hardware processor 1104B (hereinafter 
alternatively collectively referred to as processor 1104) and a 
memory 1106, such as random access memory (RAM). The 
computer 1102 may be coupled to other devices, including 
input/output (I/O) devices such as a keyboard 1114, a mouse 
65 device 1116 and a printer 1128. 

In one embodiment, the computer 1102 operates by the 
general purpose processor 1104A performing instructions 
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defined by the computer program 1110 under control of an 
operating system 1108. The computer program 1110 and/or 
the operating system 1108 may be stored in the memory 1106 
and may interface with the user and/ or other devices to accept 
input and commands and, based on such input and commands 
and the instructions defined by the computer program 1110 
and operating system 1108 to provide output and results. 

Output/results may be presented on the display 1122 or 
provided to another device for presentation or further pro- 
cessing or action. In one embodiment, the display 1122 com- 
prises a liquid crystal display (LCD) having a plurality of 
separately addressable pixels formed by liquid crystals. Each 
pixel of the display 1122 changes to an opaque or translucent 
state to form a part of the image on the display in response to 
the data or information generated by the processor 1104 from 
the application of the instructions of the computer program 
1110 and/or operating system 1108 to the input and com- 
mands. Other display 1122 types also include picture ele- 
ments that change state in order to create the image presented 
on the display 1122. The image may be provided through a 
graphical user interface (GUI) module 1118A. Although the 
GUI module 1118A is depicted as a separate module, the 
instructions performing the GUI 118B functions can be resi- 
dent or distributed in the operating system 1 1 08, the computer 
program 1110, or implemented with special purpose memory 
and processors. 

Some or all of the operations performed by the computer 
1102 according to the computer program 1110 instructions 
may be implemented in a special purpose processor 1 104B. In 
this embodiment, some or all of the computer program 1110 
instructions may be implemented via firmware instructions 
stored in a read only memory (ROM), a programmable read 
only memory (PROM) or flash memory within the special 
purpose processor 1104B or in memory 1106. The special 
purpose processor 1104B may also be hardwired through 
circuit design to perform some or all of the operations to 
implement the present invention. Further, the special purpose 
processor 1104B may be a hybrid processor, which includes 
dedicated circuitry for performing a subset of functions, and 
other circuits for performing more general functions such as 
responding to computer program instructions. In one embodi- 
ment, the special purpose processor is an application specific 
integrated circuit (ASIC). 

The computer 1102 may also implement a compiler 1112 
which allows an application program 1110 written in a pro- 
gramming language such as COBOL, C++, FORTRAN, or 
other language to be translated into processor 1104 readable 
code. After completion, the application or computer program 
1110 accesses and manipulates data accepted from I/O 
devices and stored in the memory 1106 of the computer 1102 
using the relationships and logic that was generated using the 
compiler 1112. 

The computer 1102 also optionally comprises an external 
communication device such as a modem, spacecraft link, 
Ethernet card, or other device for accepting input from and 
providing output to other computers. 

In one embodiment, instructions implementing the operat- 
ing system 1108, the computer program 1110, and/or the 
compiler 1112 are tangibly embodied in a computer-readable 
medium, e.g., data storage device 1120, which could include 
one or more fixed or removable data storage devices, such as 
a zip drive, floppy disc drive 1124, hard drive, CD-ROM 
drive, tape drive, or a flash drive. Further, the operating sys- 
tem 1108 and the computer program 1110 are comprised of 
computer program instructions which, when accessed, read 
and executed by the computer 1102, causes the computer 
1102 to perform the steps necessary to implement and/or use 
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the present invention or to load the program of instructions 
into a memory, thus creating a special purpose data structure 
causing the computer to operate as a specially programmed 
computer executing the method steps described herein. Com- 
5 puter program 1110 and/or operating instructions may also be 
tangibly embodied in memory 1106 and/or data communica- 
tions devices 1130, thereby making a computer program 
product or article of manufacture according to the invention. 
As such, the terms “article of manufacture,” “program storage 
to device” and “computer program product” or “computer read- 
able storage device” as used herein are intended to encompass 
a computer program accessible from any computer readable 
device or media. 

Of course, those skilled in the art will recognize that any 
15 combination of the above components, or any number of 
different components, peripherals, and other devices, may be 
used with the computer 1102. 
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What is claimed is: 

1. A method of computing a thrust profile to land a space- 
craft at or near a surface target, wherein the spacecraft is 30 
subject to a non-convex thrust constraint, comprising the 
steps of: 

computing a minimum landing error descent solution for 
the spacecraft subject to a convexified non-convex con- 
straint, wherein the descent solution includes a landing ^ 
error range; and 

determining the thrust profile for the spacecraft subject to 
the non-convex constraint from the computed minimum 
landing error descent solution for the spacecraft subject 
to the convexified non-convex constraint. 

2. The method of claim 1, wherein the step of determining 40 
the thrust profile for the spacecraft subject to the non-convex 
constraint from the computed minimum landing error descent 
solution subject to the convexified non-convex thrust con- 
straint comprises the steps of: 

computing a minimum fuel descent solution for the space- 45 
craft subject to the convexified non-convex constraint 
and the landing error range; and 
determining the thrust profile for the spacecraft subject to 
the non-convex constraint from the computed minimum 
fuel descent solution for the spacecraft subject to the 50 
convexified non-convex constraint. 

3. The method of claim 2, wherein the step of computing a 
minimum landing error descent solution for the spacecraft 
subject to the convexified non-convex thrust constraint com- 
prises the steps of computing 
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wherein 

r(t) is a position vector of the spacecraft as a function of 
flight time; 

r(t) is a velocity vector of the spacecraft as a function of 
flight time; 

r(t) is an acceleration vector of the spacecraft as a function 
of flight time; 

t is the time from beginning of powered landing; 
t(0) is the time at the beginning of powered landing 
tffjAyis the time at the ending of powered landing; 

T c (t) is a net thrust force vector acting on the spacecraft as 
a function of flight time; 

m(t) is a mass of the spacecraft as a function of flight time; 
rii(t) is a mass change of the spacecraft as a function of 
flight time; 

X is the set of feasible positions of the spacecraft; 
m(0) is the spacecraft mass at the beginning of the thrust 
profile; 

nifty) is the spacecraft mass at the end of the thrust profile; 
m Aj , is the initial mass of the spacecraft without fuel; 
m u , et is the initial mass of the spacecraft with fuel; 
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wherein \ sp is the specific impulse of the thruster, g e is the 
earth’s gravitational constant, Tj and T, are the lower 
and upper limits of the thrust force that can be provided 
by each thruster; 

p^nTj cos <[> is the minimum thrust value available from 
the thruster(s); 

p 2 =nT 2 cos (j) is the maximum thrust value available from 
the thrusters(s); 

e, is a unit vector of all zeros except the \ th row, which is 
unity; and 

r(t)=slack variable that bounds thrust magiitude. 

4. The method of claim 3, wherein 
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subject to the first set of constraints is computed as 
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subject to a second set of constraints comprising: 
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subject to the first set of constraints comprising the con- 
straints of: 

r(t)=g+ T c (t)/m(f) , m(t)=aT{t) 
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wherein g is the acceleration of gravity of the planet or 
object the spacecraft is landing on. 

5. An apparatus for computing a thrust profile to land a 
spacecraft at or near a surface target, wherein the spacecraft is 

5 subject to a non-convex thrust constraint, comprising the 
steps of: 

a spacecraft processor; 

a memory, communicatively coupled to the processor, the 
memory for storing instructions comprising 
to instructions for computing a minimum landing error 
descent solution for the spacecraft subject to a con- 
vexified non-convex constraint, wherein the descent 
solution includes a landing error range; and 
instructions for determining the thrust profile for the 
15 spacecraft subject to the non-convex constraint from 

the computed minimum landing error descent solu- 
tion for the spacecraft subject to the convexified non- 
convex constraint; 

a thruster, communicatively coupled to the spacecraft pro- 
20 cessor, for maneuvering the spacecraft according to the 
determined thrust profile. 

6. The apparatus of claim 5, wherein the instructions for 
determining the thrust profile for the spacecraft subject to the 
non-convex constraint from the computed minimum landing 

25 error descent solution subject to the convexified non-convex 
thrust constraint comprise: 

instructions for computing a minimum fuel descent solu- 
tion for the spacecraft subject to the convexified non- 
convex constraint and the landing error range; and 
30 instructions for determining the thrust profile for the space- 

craft subject to the non-convex constraint from the com- 
puted minimum fuel descent solution for the spacecraft 
subject to the convexified non-convex constraint. 

7. The apparatus of claim 6, wherein the instructions com- 
35 puting a minimum landing error descent solution for the 

spacecraft subject to the convexified non-convex thrust con- 
straint comprises the instructions for computing 


wherein u(t A )=is the control input, and 


o-fe) = 


i~to) 

m(h) 


p 1 (t)=p 1 e z °, wherein z 0 =a lower bound on In m(t) 
i l 2 (t) = P e ~ r °) wherein z 0 =a lower bound on In m(t) 


40 min ||r(r / )|| 2 

rf,T c o).m 


subject to the first set of constraints comprising the con- 
straints of: 

f(t)=g+T c (t)/m(t), rh(t)=aT(t) 


yt 


r(tt) 
Kh) 
z(t t ) 


||r„(f)sr(r), o< P iSr(f)<p 2 

r{t)(EX<J te.fi, t f \ 
m(0)=m we „ m(tj)^m dry 


r(t A .) is r(t)t A 

rfe) isr(t)att A . 55 

z(t A ) is z(t) at t=t A 

e 4 r =A transpose of a vector of all zeros except the fourth 
element, which is unity 
y N r A transpose of y k at the final time step N 
<J> A , A a , 'Pj., Y /f are discrete time state transition matrices 60 
describing the solution to r(t)=u(t)+g and 



r( 0)=r„, r(0)=r o 

rififerO, rit Jr ° 

wherein 

r(t) is a position vector of the spacecraft as a function of 
flight time; 

r(t) is a velocity vector of the spacecraft as a function of 
flight time; 

r(t) is an acceleration vector of the spacecraft as a function 
of flight time; 

t is the time from beginning of powered landing; 

t(0) is the time at the beginning of powered landing 

t(t)=tyis the time at the ending of powered landing; 



US 8,489,260 B2 


27 


T c (t) is a net tlirust force vector acting on the spacecraft as 
a function of flight time; 

m(t) is a mass of the spacecraft as a function of flight time; 
m(t) is a mass change of the spacecraft as a function of 
flight time; 

X is the set of feasible positions of the spacecraft; 
m(0) is the spacecraft mass at the beginning of the tlirust 
profile; 

m(h) is the spacecraft mass at the end of the tlirust profile; 
m d is the initial mass of the spacecraft without fuel; 
m wel is the initial mass of the spacecraft with fuel; 


/ sp g,COS0’ 


15 


earth’s gravitational constant, Tj and T 2 are the lower 
and upper limits of the tlirust force that can be provided 
by each thruster; 

p^nTj cos <() is the minimum thrust value available from 
the thruster(s); 

p 2 =nT 2 cos <(> is the maximum tlirust value available from 
the thrusters(s); 

e, is a unit vector of all zeros except the \ ,h row, which is 
unity; and 

r(t)=slack variable that bounds thrust magnitude. 

8. The apparatus of claim 7, wherein 
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25 
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-continued 
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wherein u(t t )=is the control input, and 
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35 
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elr k ti£it 2 (t)[l-(Fy k — Zo(r*))l 


yt = 


r(ft) 
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r (ft) is 

r(t t ) is (t) at X k 
z(y is Z(t) at t=t* 

e 4 r =A transpose of a vector of all zeros except the fourth 
element, which is unity; 
y N T =A transpose of y k at the final time step N 
(j> t . Y k are discrete time state transition matrices 

describing the solution to r(t)=u(t)+g and 


Mt) 

i= rtF) = - aa{t) 


wherein g is the acceleration of gravity of the planet or 
object the spacecraft is landing on. 

50 9. An apparatus for computing a tlirust profile to land a 

spacecraft at or near a surface target, wherein the spacecraft is 
subject to a non-convex tlirust constraint, comprising: 
means for computing a minimum landing error descent 
solution for the spacecraft subject to a convexified non- 
convex constraint, wherein the descent solution includes 
a landing error range; and 

means for determining the tlirust profile for the spacecraft 
subject to the non-convex constraint from the computed 
minimum landing error descent solution for the space- 
craft subject to the convexified non-convex constraint. 

60 10. The apparatus of claim 9, wherein the means for deter- 

mining the tlirust profile for the spacecraft subject to the 
non-convex constraint from the computed minimum landing 
error descent solution subject to the convexified non-convex 
tlirust constraint comprises: 

65 means for computing a minimum fuel descent solution for 
the spacecraft subject to the convexified non-convex 
constraint and the landing error range; and 
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means for determining the thrust profile for the spacecraft 
subject to the non-convex constraint from the computed 
minimum fuel descent solution for the spacecraft subject 
to the convexified non-convex constraint. 

11. The apparatus of claim 10, wherein the means for 
computing a minimum landing error descent solution for the 
spacecraft subject to the convexified non-convex thrust con- 
straint comprises means for computing 


min ||r(r/)|| 2 
if.Tci-mr 


30 

12. The method of claim 11, wherein 

min ||r(r / )|f 
tf, TAW)" 11 

subject to the first set of constraints is computed as 
minll^ll 2 , 

io subject to a second set of constraints comprising: 

||£V,y*77|| < ejtii; for k = 0 /V 


subject to the first set of constraints comprising the con- 
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wherein 
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r(tj) T e 1=0, r(t f )= 0 

wherein 

r(t) is a position vector of the spacecraft as a function of 
flight time; 

r(t) is a velocity vector of the spacecraft as a function of 
flight time; 

r(t) is an acceleration vector of the spacecraft as a function 
of flight time; 

t is the flight time of the spacecraft 

t(0) is the time at the beginning of powered landing 

^fH^is the time at the ending of powered landing; 

T c (t) is a net thrust force vector acting on the spacecraft as 
a function of flight time; 

m(t) is a mass of the spacecraft as a function of flight time; 
m(t) is a mass change of the spacecraft as a function of 
flight time; 

X is the set of feasible positions of the spacecraft; 
m(0) is the spacecraft mass at the beginning of the thrust 
profile; 

m(ty) is the spacecraft mass at the end of the thrust profile; 
a\ dty is the initial mass of the spacecraft without fuel; 
m MI( , r is the initial mass of the spacecraft with fuel; 
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wherein p 7 is a vector of parameters and p^GR 4 
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wherein u^j^s the control input, and 


1 
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wherein \ sp is the specific impulse of the thruster, g e is the 
earth’s gravitational constant, Tj and T 2 are the lower 
and upper limits of the thrust force that can be provided 
by each thruster; 60 

p^nTj cos <|) is the minimum thrust value available from 
the thruster(s); 

p 2 =nT, cos <|) is the maximum thrust value available from 
the thrusters(s); 

e ; is a unit vector of all zeros except the i rh row, which is 65 
unity; and 

T(t)=slack variable that bounds thrust magnitude. 
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rfe) is (t) at t k 
z(t*) is z(t) at t=H k 
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e 4 r =A transpose of a vector of all zeros except the fourth 
element, which is unity 
y N T =A transpose of y k at the final time step N 
<b k , A k , 'P t , Y k are discrete time state transition matrices 
describing the solution to r(t)=u(t)+g and 5 


z 


m(r) 

m(t) 


= -acr(t) 


10 


wherein g is the acceleration of gravity of the planet or 
object the spacecraft is landing on. 



